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ABSTRACT 

We present an extragalactic perspective of an extended stellar tidal stream 
wrapping around the edge-on, spiral galaxy NGC 5907. Our deep images reveal 
for the first time a large scale complex of arcing loops that is an excellent example 
of how a low-mass satellite accretion can produce an interwoven, rosette-like 
structure of debris dispersed in the halo of its host galaxy. The existence of 
this structure, which has probably formed and survived for several Gigayears, 
confirms that halos of spiral galaxies in the Local Universe may still contain a 
significant number of galactic fossils from their hierarchical formation. 

To examine the validity of the external accretion scenario, we present N- 
body simulations of the tidal disruption of a dwarf galaxy-like system in a disk 
galaxy plus dark halo potential that demonstrate that most of the observed tidal 
features observed in NGC 5907 can be explained by a single accretion event. 
Unfortunately, with no kinematic data and only the projected geometry of the 
stream as constraint, the parameters of our model are considerably degenerate 
and, for now, must be considered illustrative only. 

Interestingly, NGC 5907 has long been considered a prototypical example of 
a warped spiral in relative isolation. The presence of an extended tidal stream 
challenges this picture and suggests that the gravitational perturbations induced 
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by the stream progenitor must be considered as a possible cause for the warp. 
The detection of an old, complex tidal stream in a nearby galaxy with rather 
modest instrumentation points to the viability of surveys to find extragalactic 
tidal substructures around spiral galaxies in the Local Volume (< 15 Mpc) - 
with the prospect of obtaining a census with enough statistical significance to be 
compared with cosmological simulations. 

Subject headings: galaxies: individual (NGC5907) — galaxies: dwarf — galaxies: 
evolution — galaxies: interactions — galaxies: halos — dark matter 

1. Introduction 

The Cold Dark Matter (ACDM) paradigm predicts that stellar halos surrounding large 
spiral galaxies formed through the accretion and tidal disruption of satellite galaxies, a notion 
previously postulated (Searle & Zinn 1978) on empirical grounds from the character of stellar 
populations found in our own Milky Way halo. Simulations have shown that the fossil records 
of those merger events may be detected nowadays in the form of coherent stellar structures 
in the outer regions of the host galaxies. The most spectacular of these fossil structures 
are probably the formation of dynamically cold stellar streams, which are similar to long 
rivers of stars (and likely dark matter; Peharrubia et al. (2008b)) tidally stripped from a 
disrupting dwarf galaxy and wrapped around the host galaxy disk, roughly tracing the orbit 
of the progenitor satellite. The now well-studied Sagittarius (Sgr) tidal stream surrounding 
the Milky Way (Ibata et al. 2001b; Majewski et al. 2003; Martmez-Delgado et al. 2004) 
and the giant stream in Andromeda galaxy (Ibata et al. 2001a) are archetypes in the Local 
Group of these kinds of fossilized satellite galaxy mergers and provide obvious support for 
the scenario that tidally disrupted dwarf galaxies are an important contributor to the stellar 
halos of disk galaxies. 

Unfortunately, while the Milky Way and even M31 streams are close enough to be 
resolved into individual stars that can be studied in detail, these nearby streams can also span 
vast areas of sky, which currently presents particular challenges to their study. In contrast, 
tidal streams around more distant systems can be imaged in commensurately smaller fields 
of view, and — at least for those streams seen in diffuse light — with equal contrast, given 
the constancy of surface brightness with distance. However, over the past decade, only a few 
cases of confirmed stellar tidal streams have been detected in disk galaxies outside the Local 
Group (e.g., Malin & Hadley 1997; Forbes et al. 2003; Pohlen et al. 2004, and references 
therein), in contrast with the increasing number of tidal structures detected in the Milky 
Way and M31 in recent years. The presence of a multitude of tidal streams, arcs, shells and 
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other irregular structures possibly related to ancient merger events in the deep panoramic 
view of the Andromeda halo (Ibata et al. 2007) reveals the level of stellar sub-structure that 
might be present in the halos of nearby external spiral galaxies. 

Extragalactic streams also offer a unique perspective to address some interesting ques- 
tions: Is the abundant number of stellar streams exceptional in the Local Group or are tidal 
streams as common as has been predicted by current cosmological models? What can we 
learn about the dark-matter halos of the host galaxies from these streams? Can we use these 
systems to study the star formation history of the disrupting dwarf galaxies? The frequency 
of streams, their stellar population and orbital properties provide important clues to the na- 
ture of the stream progenitors, lend insights into the underlying gravitational potential and 
structure of the massive dark halo they inhabit, and constrain hierarchical models of galaxy 
formation favored by current cosmological models. In that respect, an interesting parameter 
that may possibly be constrained through the detection of tidal streams in external galaxies 
is the flattening of the dark matter halos surrounding the host galaxies. The spatial and 
kinematic distributions of stream particles are fairly sensitive to this quantity (e.g. Johnston, 
Choi, & Guhathakurta 2002; Helmi 2004; Johnston et al. 2005; Law, Johnston, & Majewski 
2005; Peharrubia et al. 2006, and references therein) and the shapes and orientations of dark 
matter halos with respect to the disk axi-symmetry plane reflect the nature of dark matter 
(see, e.g. Dubinsky et al. 1991; Yoshida et al. 2000; Dave et al. 2001, for shape estimates 
for cold, self-interacting and hot dark matter models, respectively). For example, ACDM 
models predict that dark matter halos are triaxial, with minor axes perpendicular to the 
disk plane. On average the minor-to- major axis ratio near the center is c/a ~ 0.7 (Hayashi, 
Navarro, & Springel 2007), although halos may become considerably rounder after the bary- 
onic galactic components have formed (Kazantzidis et al. 2004) or if dark matter is warm 
(Biermann & Munyaneza 2007). Thus, using tidal streams in external galaxies might be an 
ideal tool to obtain a significant statistical sample of flattening values that can be compared 
against predictions from cosmological models. 

These arguments suggest that a systematic survey for tidal streams around extragalactic 
systems would provide a new way to explore galaxy structure and evolution, with direct 
constraints on cosmological models. Undertaking such a tidal stream survey was urged by 
Johnston, Sackett, & Bullock (2001) in their comprehensive discussion of the appearance 
and detectability of extragalactic tidal streams. Their analysis suggests that a survey of 100 
parent galaxies reaching to a surface brightness of 30 mag arcsec -2 would reveal many tens 
of tidal features; a slightly higher estimate of about one per galaxy is obtained in their most 
recent reanalysis of the problem (Bullock & Johnston 2005). These authors also point out 
the advantage of exploring thin, edge-on disk systems because of the favorable orientation for 
detecting low surface brightness features against light of the parent system and the ability 
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to assess limits on the mass ratio of the satellite and parent galaxy by whether the disk of 
the parent appears to be disturbed. 

Johnston, Sackett, & Bullock (2001) estimated that a search strategy designed to obtain 
reliable photometry along the stream length would require about ten times longer exposure 
times that those aimed at detection only. We adopted this approach of striving for deep 
images to detect the full path of the tidal stream, which, in the absence of kinematic data, is 
the main input for comparison with N-body simulations (see Sec. 4). Using this strategy, we 
started a pilot survey on some nearby galaxies which has proven the concept with the dis- 
covery of a giant stellar tidal stream around the spiral galaxy NGC4013 (Martinez- Delgado 
et al. 2008) and the detection of the faintest parts of the extended disk of the nearby spiral 
galaxy M94 (Trujillo et al. 2008, in preparation). 

With the hope of improving our understanding of the known tidal stream in NGC 5907, 
discovered by Shang et al. (1998), and to demonstrate the sensitivity of the small aperture 
telescope used in our survey (see Sec. 2. 2), we reobserved this system as a commissioning 
target. In this paper we present our very deep image of NGC 5907 that reveals additional 
loops of debris in the halo (see Sec. 3.1) beyond those reported by Shang et al. (1998). N-body 
simulations, described in Sec. 4 show that most of the newly found features can be explained 
as pieces of a spectacular, multiply-wrapped stream of debris from the tidal disruption of 
a single companion galaxy. We also discuss in Sec. 5 the applications of such surveys to 
glean a deeper theoretical understanding of galaxy formation and evolution, beyond placing 
strong constraints on present cosmological scenarios. This foray into a more systematic look 
for tidal streams in the nearby universe has yielded promising results that encourage more 
aggressive attention to this new way of understanding galaxy structure and evolution as 
advocated by the above theoretical studies. 

2. Target, Observations and Data Reduction 

2.1. The Target: NGC 5907 

NGC 5907 is a nearby (14 Mpc; Zepf et al. 2000), edge-on, Sc galaxy, very similar 
(e.g., type, absolute magnitude, and rotation speed) to the Milky Way. NGC 5907 has 
frequently been the target of deep surveys (radio, optical and near-infrared) because it is 
an apparently undisturbed, almost edge-on disk useful for exploring the properties of spiral 
disks, determining the properties of dust and gas in disk galaxies, testing for the prevalence 
of luminous thick disks and halos, searching for warps, and exploring the strength of dark 
matter in disks (e.g. Sancisi 1976; van der Kruit & Searle 1981; Casertano 1983; Skrutskie, 
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Shure, & Beckwith 1985; Sasaki 1987; Barnaby & Thronson 1992; Morrison, Boroson, & 
Harding 1994; Sackett et al. 1994; Lequeux 1996, 1998; Rand 1996; Dumke et al. 1997; Rudy 
et al. 1997; James & Casali 1998; Xilouris et al. 1999; Zheng et al. 1999). Although VLA HI 
maps of this galaxy showed that it has a large-scale warped gas disk (Sancisi 1976), NGC 5907 
was long considered to be a prototypical example of a warped disk galaxy existing in relative 
isolation since the identification of a neighbouring, interacting galaxy was lacking. However, 
deep imaging of NGC 5907 by Shang et al. (1998) has revealed a very faint, elliptically- 
shaped ring around the disk, and this has been interpreted to be the remnant of a tidally 
disrupted dwarf satellite galaxy, an analog to the debris stream created by the Milky Way-Sgr 
interaction, although at higher surface brightness. 



2.2. Observation and Data Reduction 

We have obtained deep optical images of NGC 5907 with the 0.5-meter Ritchey- Chretien 
telescope of the BlackBird Remote Observatory (BBRO) situated in the Sacramento Moun- 
tains (New Mexico, USA). We used a Santa Barbara Instrument Group (SBIG) STL-11100 
CCD camera, which yields a large field of view (27.7' x 18.2') at a plate scale of 0.45 arcsec- 
onds pixel -1 . The image set consists of multiple deep exposures with a non-infrared, clear 
luminance (3500 < A < 8500), a red, a green and a blue filter from the SBIG Custom Scien- 
tific filter set. Table 1 provides a summary of the images collected during different dark sky 
observing runs during the period June through August 2006. Column 3 refers to the total 
exposure time of the co-added images from each filter obtained with each run. 

The science images were reduced using the standard procedures for bias correction and 
flat-fielding. A master dark and bias frame was created by combining 10 dark sub-exposures 
each produced at the same exposure length and camera temperature settings used for the 
luminance and the filtered images. A master flat was produced by combining 10 separate 
sky flat exposures for each filter. The red, green and blue filtered exposures were separately 
combined (using a median procedure) to produce red, green and blue master images that 
represented the total exposure time for each color filter. These color images were obtained 
to investigate the color of the structures surrounding the galaxy and to extend the overall 
exposure of the final luminance image. The master images for each filter were subsequently 
summed to produce a synthetic luminance image that represented the total exposure time of 
all the color filter data. The synthetic luminance image was then combined (using a median 
procedure) with all the clear filtered luminance sub-exposures (see Table 1) to increase the 
contrast of the diffuse light. The resulting final image thus represents the sum of all available 
CCD exposures collected for this project, with an accumulated exposure time of 11.35 hours 
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( including 5.75 hours in the clear luminance filter). 

To enhance the faint structures around NGC 5907 and disentangle their true appearance 
from image noise effects, the contrast and detail of the summed image was optimized using 
two well-known image processing techniques: i)a histogram equalization of the image using a 
non-linear transfer function 1 ; and ii) Gaussian blur filtering (Davies 1990; Haralick & Shapiro 
1992). The first method is very effective to reveal very faint features in images that contains 
a large dynamical range, as it is the case of a spiral galaxy with a very bright, compact 
stellar disk surrounded by a very faint tidal stream. With the purpose of suppressing the 
bright portions of the image and intensifying the fainter parts of the stream, histogram 
equalization employs a monotonic, non-linear mapping which re-assigns the intensity values 
of pixels in the input image such that the output image contains an uniform distribution of 
intensities. In our case, we used an iterative process that involves several passes of a S-shaped 
histogram operator to the total image. While this method was successful in revealing the 
presence of the surrounding faint structure, it also resulted in the over saturation of contrast 
and significantly reduced intensity variances within the stellar stream. Once the image was 
contrasted, a low-pass Gaussian filter with a radius of 2 pixels was applied to reduce the low 
level noise contaminating the image. The resulting noise-filtered image is shown in Fig.l. 

As discussed above, our search technique was designed to very clearly reveal the position 
and morphology of faint structures (for example, by combining all the available images 
obtained in different filters with the clear luminance deep images), but does not permit 
accurate photometry. Therefore, follow-up photometric observations are needed to measure 
accurate surface brightnesses and to estimate the total luminosity of the detected tidal debris 
(see Sec. 3.3). Furthermore, the use of an uncalibrated luminance filter prevents us from 
estimating the surface brightness limit that we reached in our deep probe of the NGC 5907 
halo in any particular, standard photometric band. However, the detection of obvious diffuse 
light structures fainter than those reported in previous studies (with a surface brightness as 
faint as Hr ~ 27.7 mag arcsec -2 ; see Sect. 3.3), suggests that our images reach a surface 
brightness magnitude fainter than the R = 27 mag arcsec -2 limit of the deep NGC 5907 
image by Morrison, Boroson, & Harding (1994) and Sackett et al. (1994), and the R = 28.7 
mag arcsec -2 limited image by Zheng et al. (1999). 

Regardless, this result is a striking example of the scientific potential of modest aperture 
telescopes (0.5- meter) operating under very dark skies for probing the diffuse structure of 
galaxies. 



1 see http://http://homepages. inf.ed.ac.uk/rbf/HIPR2/histeq. htm and references within 
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3. The Tidal Stream(s) of NGC 5907 
3.1. Observed Morphology of the NGC 5907 System 

3.1.1. Truncation and Optical Warp 

In this Section we focus on the outer, diffuse light properties revealed in Figures 1 and 
2, but several comments about the disk are worth pointing out briefly. 

Our image shows that the disk of NGC 5907 extends in diameter well beyond where 
most previous images have tracked its disk light. The existence of this light challenges 
previous claims that the disk is sharply truncated (in the sense of the galaxy's outer edge) 
at a nominal cut-off radius of about 6' from the center (van der Kruit & Searle 1981; 
Sasaki 1987; Morrison, Boroson, & Harding 1994). From the original radial profiles (van der 
Kruit & Searle 1981; Sasaki 1987, Fig.9 and Fig.7, respectively) there is clearly a radial 
break feature visible, however, it should be better described as a sharp change of slope than 
a complete cut-off. This profile shape, well fitted with a broken exponential function, is a 
relatively common feature of late-type spiral galaxies (Pohlen & Trujillo 2006; Erwin, Pohlen, 
& Beckman 2008). 

Figures 1 and 2 clearly show that NGC 5907 has a stellar warp, which clearly extends 
after the nominal cut-off radius. The observable faint radial extensions of the disk are 
strongly bent in the same sense as the obvious HI warp found by Sancisi (1976), and are 
most strongly visible on the northwest disk edge. Note, these same extended disk features 
can also be seen in Fig. 2 of Shang et al. (1998), and just about in the near infrared image 
of Barnaby & Thronson (1992). Evidence for the warp in optical light was also previously 
seen in the optical disk at smaller radii by Sasaki (1987), Florido et al. (1992), and Morrison, 
Boroson, & Harding (1994) and was, of course, pointed out as well by Shang et al. (1998). 

3.1.2. Diffuse Light Structures 

The most striking features, and the primary focus of the present discussion, is the newly 
revealed network of arcing loops of diffuse light around NGC 5907 that are likely the product 
of the classical satellite tidal disruption process. Our deep image (see Fig. 2) reveals that 
the ring-shaped feature to the southeast, we refer to this here as the SE loop, discovered 
by Shang et al. (1998) is only the brighter part of a large-scale complex of stellar debris 
arcs around NGC 5907. We count about a half-dozen other additional arc-like features in 
our image. On the eastern side, it is possible to detect the full extension of the SE loop, 
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turning around and falling toward the disk. This yields a vertical structure (labeled El in 
Fig. 2) that crosses the galactic plane and possibly extends to the other side of the galaxy. 
In addition, an independent structure (labeled E2), somewhat parallel to the SE-loop, can 
be seen to be shifted about ~ 10 ' to the northwest of the Shang et al. (1998) structure and 
sweeping almost 180° around the center of NGC 5907. 

Although the accretion origin of this structure appears a reasonable hypothesis, we must 
point out its extraordinary complexity — especially on the western side of NGC 5907 where 
we find a faint, long tail extending ~ 25 kpc from the disk, labeled SW tail in Fig. 2, and 
previously reported only in some studies (see Fig. lb in Reshetnikov & Sotnikova 2000). 
From south to north, three additional arcs are detected on the west side of the disk (labeled 
Wl, W2 and W3, respectively). These features appear to cross each other fairly closely to 
the disk major axis of NGC 5907, which makes it difficult to discern which of these arcs are 
the continuation of the features observed on the eastern side, that is, the SE loop, El, and 
E2. 

For example, it is obvious that the SE loop extends to the western side, but it becomes 
unclear which path it takes (e.g., SW tail or Wl). On the other hand, it is clear that the W3 
arc is the westernmost continuation of the E2 loop. Interestingly, this stream piece turns 
around, crosses the W2 arc and finishes in the main disk again, forming the Wl arc. It is 
however more difficult to assess whether the W2 arc (the most diffuse discernible structure, 
highly contaminated by the diffuse light of the galaxy disk) or Wl could be the western 
extension of the El tail. 



3.1.3. Evidence for a Single Coherent Structure 

With the image in Fig. 2 alone it is impossible to conclude whether the smaller features 
on the western side are the continuation of the loops on the eastern side and so whether 
the whole set of these identified tidal structures originate from a single or different accretion 
events. In general, the largest features (e.g., the SE and E2 loops) have similar shapes, 
angular lengths, radius of curvature, and separations from the center of NGC 5907. In 
addition, all of the features (large and small) share a similar projected width, ~ 7 kpc 
(consistent with the value reported by Shang et al. 1998, see Sec. 3. 3), as well as a similar, 
rather bright surface brightness (see Sec. 3. 3). These similarities all provide circumstantial 
evidence that the various features may all come from a single, old 2 stellar tidal stream - 



2 Following Johnston, Sackett, & Bullock (2001), the age of the trail here refers to the epoch when the 
debris were unbound from the satellite core (see Fig. 4 ) and not to the age of its stellar population, that is 
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similar to the multiply-wrapped Sgr dwarf galaxy stream around the Milky Way. These 
structures provide a rather unique and spectacular external view of the classical rosette-like, 
multi-loop patterns expected for tidally disrupting dwarf galaxies seen in N-body simulations 
(such as those made of the Sgr system — e.g., Johnston, Spergel, & Hernquist 1995; Ibata 
& Lewis 1998; Martinez-Delgado et al. 2004; Law, Johnston, & Majewski 2005). There are 
further analogies of the NGC 5907 structures with the Sgr stream in that they have a similar 
linear cross-sectional width (suggesting similarly massed progenitors), apogalacticon and 
apparent orbital ellipticity (suggesting a similar orbit), and surface brightness and length 
(suggesting a similar mass loss rate and age of interaction). In Section 4 we explore this 
working hypothesis in detail and provide a comprehensive, illustrative model that further 
supports the notion that a single, Sgr-like system could create a substantial fraction of the 
newly found diffuse structures. 

Interestingly, while the surface brightness of these loops does show some modulation, 
nowhere among these features do we see a bright spot that could be identified with a re- 
maining dwarf satellite core. Thus, our image provides no definitive insights on the current 
position and final fate (that is, whether or not it is completely disrupted by now) of the 
progenitor satellite galaxy (since it could be hidden behind the disk). However, the full 
panoramic view of the NGC 5907 system in Figure 1 provides no evidence for any associa- 
tion of the stellar streams with the dwarf irregular galaxy PGC 54419 (a confirmed satellite 
of NGC 5907 situated at 36.9 kpc projected distance from its center), as proposed by Shang 
et al. (1998). We shall return to this subject in Sec. 5. 



3.2. Diffuse Light: Halo or Stream? 

As mentioned above, there have been a number of deep photometric studies of NGC 5907, 
and a popular point of discussion is the claimed detection and properties of " excess" , non- 
disk light away from and perpendicular to the disk and commonly attributed to a diffuse 
luminous halo (Sackett et al. 1994; Lequeux 1996, 1998; Rudy et al. 1997; James & Casali 
1998; Zepf et al. 2000; Irwin & Madden 2006). Figure 1 clearly shows how these previous 
studies may have been affected by the coherent tidal stream structures, especially Wl and 
W2, which lie along the southwestern minor axis at distances where such studies have made 
the "halo" detections. Indeed, the very shape of the "inner triangle" of light formed by Wl 
and W2 can be discerned in the image analyzed by Lequeux (1996, see, e.g., their Fig. 1), 
and a similar asymmetry seen in the photometric data by Sackett et al. (1994) may also be 



unknown. 
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explained by these same structures (e.g., see the lower right panel of Fig. 2 in Sackett et al. 
1994). 

One wonders if the presence of these structures may also explain the differences in the 
colors of the "halo" light found on either side of the disk shown in Figure 1 of Lequeux 
(1998). Zheng et al. (1999) have previously suspected that the presence of tidal debris 
light (particularly on the western side) might be significantly contaminating the "halo" light 
that was previously reported; these authors concluded that, when combined with residual 
scattered light from foreground stars, this debris effectively reduces the need for a faint, 
extended halo. We concur that indeed these previous studies were probably strongly affected 
by the presence of the "W" features — but of course this is actually "halo" light, it is just 
in the form of tidal debris in what is obviously the active hierarchical formation of the 
NGC5907 halo. However, a low surface brightness halo (Er ~ 30-31 mag arcsec -2 ) similar 
to M31 one could be present (Fliri et al. 2008, in preparation). 

3.3. Estimating the Surface Brightness, Total Luminosity and Stellar Mass of 

the NGC 5907 Tidal Stream 

Ideally, one would like to use the BBRO image to measure the surface brightness of the 
observed stream. However, there are some concerns that prevent us from extracting this 
information directly from this image. First of all, the image is the result of the combination 
of multiple images in different filters with those taken using an extremely wide (luminance) 
filter to maximize the photon collection. This makes a calibration of the stars present in the 
BBRO image using the same stars observed through standard filters non-trivial. Significant 
color corrections are expected. 

In addition, the resulting background on the BBRO image is slightly patchy (see Fig. 1), 
and, more importantly, the field of view is not large enough to allow us to properly estimate 
the sky level. For the above reasons we have decided to adopt the value reported by Zheng 
et al. (1999) for the mean surface brightness of the stream in the NE loop in the i?-band - 
i.e. 27.7 mag arcsec -2 — and rely on the BBRO image only to determine the morphology, 
position and width of the stream in the sky. 

In order to determine the total luminosity of the tidal stream we first calculate the total 
area covered by the tidal tails. To do so, we estimate the stream width in two different strips 
A and B located in the bright stream piece detected by Sheng et al. (1998) (see Fig. 2). 
The surface brightness distribution in these two strips is shown in Fig. 3. The width of the 
stream in each of these strips is estimated using a Gaussian function and the Full Width Half 
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Maximum (FWHM) of these Gaussian fits are used as a measurement of the stream width. 
To be conservative, we take the minimum of the FWHM measured in both regions (i.e. ~103 
arcsec; see Fig. 3) as representative of the width of the stream. At the distance of NGC 5907 
the above value corresponds to ~7 kpc. Subsequently, we calculate the length of the stream 
by placing elliptical annuli at the inner at outer sides of each loop. Here we consider that 
the stream of Fig. 2 can be described by the loop LI formed by the SW-SE-E1-W2 arcs 
and the half-loop L2 formed by the W3-E2 arcs. The area of each loop is then calculated 
as A = 7rg(a 2 — ft 2 ), where a and ft are the semi-major axis of the elliptical annuli and are 
related as 6 = a — FWHM. The minor-to-major axis ratio of the ellipses has an estimated 
value of q ~ 0.5 for LI and L2. We find that for LI a = 540 arcsec= 36.6 kpc, so that 
Ali — 728 kpc 2 , and for L2 a' = 27.5 kpc, which implies Al2 = 527 kpc 2 . On doing this, the 
approximate area subtended by these pieces of the stream is A = A L1 + A L2 /2 ~ 990 kpc 2 . 

Using this estimate of the area of the SE loop and E2 loop, the reported mean surface 
brightness of 27.7 mag arcsec -2 by Zheng et al. (1999) (i.e. ~0.217 L pc -2 in the _R-band) 
implies a total i?-band luminosity for the stream of ~2xl0 8 L Q . This value is one order 
of magnitude brighter than the luminosity reported for the Sgr dwarf galaxy (~1.9xl0 7 , 
Majewski et al. 2003), but it is still consistent with the possibility that this complex structure 
of debris resulted from the disruption of a single dwarf galaxy during a recent interaction 
with NGC 5907 (see Sec. 4.2). We have also roughly estimated the total stellar mass of the 
stream by using the color R — I =0.5±0.3 given in Zheng et al. (1999). Using the prescription 
given by Bell & de Jong (2001) the above color implies a (M/L) R ~ 1.6^ 5 4 and consequently 
a total stellar mass of ~3.5xl0 8 M . This value is very uncertain due to the large error in 
the stream color. So, the total stellar mass should only be taken as a rough estimate of the 
stellar mass in the stream within an order of magnitude. Given the large uncertainties of our 
estimate it is encouraging that the above stellar mass estimate is similar to the one estimated 
(~2xl0 8 M ) by Johnston, Sackett, & Bullock (2001) using dynamical considerations. 

Unfortunately, none of the individual color filtered data set were of sufficient exposure 
length to reveal hints of the faint surrounding structures with enough signal-to-noise to 
permit detection or surface brightness measurement. Therefore, it is not possible to constrain 
a stellar population gradients along the stream from this study. 

4. N-body Simulations 

The clear picture of a repeatedly wrapped tidal stream around NGC 5907 naturally 
suggests fitting the system with N-body models. Typically the number of free parameters 
implemented by N-body models of satellite accretion in external galaxies is vast. In the 
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case of NGC5907, although the mass within the luminous radius is well constrained by 
measurements of its circular velocity curve, no information is available on the flattening of 
its potential on the scales where the tidal remnants are observed (~ 10-50 kpc in projection). 
Additionally, the present position and velocity vectors of the progenitor system are unknown 
as well as its mass and luminous profiles. 

Currently, the only available information on the stream is its surface brightness and 
projected geometry with respect to the host. This introduces a new free parameter, which 
is the projection angle between the line of sight and the orbital plane of the stream progen- 
itor. Fortunately, NGC 5907 presents an edge-on disk projection that eliminates a second 
projection angle between the line-of-sight vector and the axi-symmetry plane of the disk. Ki- 
netically data of bright dynamical tracers (as planetary nebula and globular clusters) along 
the stream would help to disentangle the projection effects and to constrain the satellite or- 
bit. Unfortunately, these cannot be readily obtained from individual stars with the present 
observational capabilities owing to the large distance of NGC 5907. 

The huge degeneracy of the problem weakens any possible constraint obtained through 
N-body modeling of the stream and, in fact, it is presently unclear how much information can 
be gleaned from streams in external galaxies using this technique. That is an open question 
that goes beyond the scope of this paper and that will be studied in a future contribution 
(Penarrubia et al., in prep). 

Nevertheless, in this paper our N-body models are only used to answer a more simple 
question: can the complex stellar structure be the result of a single merger event? This 
reduces the problem of finding one of the many possible free parameter combinations that 
yields a reasonably good match to the observations. 

4.1. Model Parameters 

We examine several illustrative N-body simulations of dwarf spheroidal systems disrupt- 
ing in a spiral galaxy having properties similar to NGC 5907. This galaxy is an Sc spiral with 
a maximum rotational velocity of V^ max = 227 km s _1 and absolute magnitude Mb — —21 
(Casertano 1983) — similar, therefore, to the Milky Way. 

Here we assume a static host galaxy potential during the satellite orbit. As Penarrubia 
et al. (2006) show, the geometry and kinematics of a tidal stream are independent of the 
past evolution of the host galaxy potential and solely reflect the present properties of the 
host. Our satellite galaxy realizations are injected in a host galaxy that has the following 
components: a Miyamoto & Nagai (1975) disk, a Hernquist (1990) bulge and a Navarro, Frenk 
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& White (1996,1997) dark matter halo. The determination of the disk, bulge and dark 
matter halo mass from the circular velocity curve of NGC 5907 is very uncertain due to the 
well-known disk-halo degeneracy. Here we adopt values that are similar to those derived 
by Just et al. (2006) by fitting the observed rotation curve: the disk and bulge masses are 
Md = 8.4 x 1O 1O M and Mb = 2.3 x 1O 1O M , respectively. The radial and vertical exponential 
scale-lengths of the disk are Rd = 6.4 kpc and Zd = 0.26 kpc, respectively. The bulge core 
radius is c = 0.6 kpc. The halo density profile can be written as 

p NFW = 5cpc (1) 
r/r s (l + r/r s ) 2 ' 

where p c =277.72 /i 2 M kpc 2 is the present-day critical density, h = 0.71 is the Hubble con- 
stant in units of 100 km s _1 Mpc -1 consistent with constraints from cosmic microwave 
background (CMB) measurements and galaxy clustering (see Spergel et al 2007 and refer- 
ences therein), 8 C is a dimensionless parameter and r s is the halo scale radius. For 8 C = 200 
we use M 2 oo = 6.8 x 1O 11 M , r s ~ 8.2 kpc and r 2 oo = 180 kpc. We must remark, however, 
that different combinations of the host galaxy parameters can yield a similar value of V Cjiaax 
(see Just et al. 2006 for a detailed disccusion). Unfortunately, this increases the degeneracy 
inherent to the models of the tidal stream surrounding NGC 5907. 

The satellite galaxy is modeled as a King (1966) profile for simplicity. More realistic, 
multi-component models that distinguish between baryonic and dark matter, such as those 
motivated by ACDM (e.g. Penarrubia et al. 2008a,b), would further increase the number 
of free parameters. Our satellite galaxies have typical parameters for nearby dwarf galaxies 
(Mateo 1998): a total mass of 2 x 1O 8 M , core radius of R c = 390 pc and King concentration 
Cfc = log 10 ( Rt/R c ) ~ 0.84, where R t is the King tidal radius. These values were selected to 
reproduce the projected width and extent of the stellar tidal stream detected in NGC 5907 
(see below). 

In order to find an illustrative model that provides a reasonable match to the projected 
distribution of debris, we run a set of N-body simulations that explore some of the satellite's 
orbital and structural parameters. Following Reshetnikov & Sotnikova (2000), we assume 
that the orbital apocenter corresponds to the maximum separation of the stream from the 
center of NGC 5907 (i.e. that of the SE loop), so that r a = 51 kpc. To estimate the orbital 
eccentricity, we vary the peri to apocentric distance ratio (r p /r a = 0.1, 0.2, 1.0). Following 
the estimates of Johnston, Sackett & Bullock (2001), we explore three different progenitor 
masses (10 8 ,2 x 10 8 and 5 x 1O 8 M ). To derive the age of the stream we integrate each 
N-body realization a maximum time of 9 Gyr, recording the positions and velocities of the 
N-body particles in 70 snap-shots (i.e. with a time interval of 0.14 Gyr). Each snap-shot is 
then visually inspected varying the projection angle (£). 
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4.2. Comparison of N-body Models to the Data 

In Fig. 4 we show an illustrative single-accretion model that produces a good overall 
match to most of the observed structures. 3 This suggests that most of the tidal loops 
could have formed from the accretion of a single dwarf satellite with an initial total mass 
of 2 x 1O 8 M . For this particular model, the satellite has lost 60% of its mass due to tidal 
stripping after 3.6 Gyr of evolution. The progenitor moves on a highly inclined orbit (orbital 
inclination i = 80° with respect to the disk plane), with relatively low eccentricity (the peri 
and apocenters are, respectively, r p = 21 kpc and r a = 51 kpc, which translates into an 
intermediate orbital eccentricity of e ~ 0.42) and a radial orbital period of T r ~ 0.9 Gyr. 
Remarkably, these parameters are fairly similar to those derived for the Sgr stream in the 
Milky Way (Law et al. 2005). In addition, we find that the projection angle formed by the 
line-of-sight and the orbital plane vector has a present value of £ ~ 57° (an edge-on view of 
the stream implies £ = 90°), which explains the considerable elongation of the stream in the 
vertical direction. 

In general, our model agrees with the previous model of Reshetnikov & Sotnikova (2000), 
although a few differences must be noticed. For example, Reshetnikov & Sotnikova (2000) 
assume a perfectly polar orbit (i = 90°), which is therefore not affected by orbital precession, 
and derive a slightly higher orbital eccentricity 0.6), probably a result of the smaller 
amount of observational constraints (only the SE loop was known at that time). 

Our model also provides a comprehensive explanation of the relation between the dif- 
ferent loops and arcs discussed in § 3.1. In particular, this solution suggests that the leading 
arm forms the SE-loop, bending toward the disk along the El loop and giving rise to the arc 
W2 after it crosses the NGC 5907 disk. The trailing tail of the tidal stream would be visible 
as the arc Wl in the west, and its continuation bends toward the east direction to form the 
W3 arc and the E2 loop after it crosses the disk. The longest stream piece that fades away 
towards the west (the SW tail; Fig. 2) would correspond to the oldest piece of the trailing 
tail with particles that, according to our model, were stripped around 3.6 Gyr ago. Although 
the shape of our illustrative stream model can reproduce most of the prominent observed 
features, it also predicts the continuation of the E2 loop that extends between the SE loop 
and El arc. This feature is not obvious in Fig.l, although the presence of an excess of light 
in the area where this arc would be located (denoted as OV in Fig. 2) is noticeable. 

Interestingly, the model also successfully reproduces the strong east-west stream asym- 



3 This figure is the final snapshot of a movie included as supplementary material in the electronic version 
of this paper. 
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metry (the SE loop reaches a projected distance of ~ 51 kpc, whereas the SW tail only 
extends out to ~ 27 kpc), which results from the eccentric orbit of the progenitor system 
and an inclined line-of-sight projection. 

This illustrative model also provides a reasonable explanation for the fate of the pro- 
genitor satellite. It suggests that the dwarf could have survived but that its remnant core 
(black points in Fig. 4) might presently be hidden behind the galaxy disk, ~ 15 kpc from the 
galactic center. If this hypothesis is correct, it leads us to conclude that the NE ring feature 
of Shang et al. (1998) corresponds to a young piece of the stream. In addition, our model 
suggests that the NGC 5907 stellar stream may be relatively old, with the fainter, outer loop 
material (SW tail in Fig. 2) becoming unbound at least 3.6 Gyrs ago (blue points in Fig. 4). 
The new tidal tails presented in this paper are considerably older than the previous stream 
pieces discovered by Shang et al. (Reshetnikov & Sotnikova 2000 estimate age ^1.5 Gyr) 
and than any definitively established piece of the Sgr tidal stream (< 2.5 Gyrs; Majewski 
et al. 2003; Law, Johnston, & Majewski 2005). The stellar stream around NGC 5907 may 
represent therefore one of the most ancient wraps of a tidal stream ever reported in the halo 
of a spiral galaxy. 



5. Discussion 

In this paper we present in detail an extragalactic perspective of a stellar tidal stream 
surrounding the nearby spiral galaxy NGC 5907. This ghostly structure is an elegant example 
of how a single low-mass satellite accretion occurring in the current epoch can produce a 
very complex, rosette-like structure of debris dispersed in the halo of its host galaxy that 
may survive for several Gigayears, a scenario predicted by N-body simulations of tidally 
disrupting stellar systems around the Milky Way (e.g., Law, Johnston, & Majewski 2005; 
Penarrubia et al. 2005; Martin et al. 2005; Martinez-Delgado et al. 2007). 

The external view of the stream wrapping NGC 5907 also offers for the first time a 
direct observational evidence that the material stripped from disrupting satellite galaxies 
(stream stars likely accompanied by dark matter) may cross the inner regions of galactic 
disks. This process will likely induce gravitational effects (e.g., creation of spiral arms, star 
formation triggering, etc) that are not yet fully understood, a phenomena that has been also 
suggested for the Sgr stream in our Galaxy (Majewski et al. 2003; Martinez-Delgado et al. 
2007) but that lacks confirmation due to our inner perspective. In addition, the presence of 
tidal streams wrapping galaxies with warped disks (see also Martinez-Delgado et al. 2008) 
may suggest satellite galaxy perturbations as the origin of those features (e.g. Velazquez & 
White 1999, Weinberg & Blitz 2006). 
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Recent cosmologically-motivated simulations have shown that in a ACDM paradigm 
the accretion of several massive (0.2-0.6M disk ) sub-halos onto the central regions of the host 
galaxy should be common from z ~ 1 to the present (Kazantzidis et al. 2007). Such accretion 
events lead to strong warping, flaring and thickening of an initially cold disk as well as to the 
formation of long-lived, ring-like stellar features in the outskirts of the disk that may locate 
several kiloparsecs off the disk plane and have surface brightnesses in the range of 25-30 
mag arcsec -2 (see Fig. 6 of Kazantzidis et al. 2007). Given the similarities in shape and 
mass between NGC 5907 and the Milky Way, it is interesting to note that our deep images 
show none of the predicted features besides a classical warped disk. A statistically relevant 
sample of deep photometric images of nearby, edge-on spiral galaxies might reveal whether 
the absence of strongly perturbed old disks poses a problem for our theoretical understanding 
of how disk galaxies form in a ACDM context. 

Our survey reveals that the presence of complex structures from ancient streams in the 
outskirts of spiral galaxies in the Local Volume ( ^ 15 Mpc) must be taken into account to 
correctly interpret the results of deep, pencil-beam photometric probes designed to constrain 
the formation and composition of galactic halos. This concern is justified in the case of 
NGC 5907, where several authors (e.g., Sackett et al. 1994; Lequeux 1996) have reported the 
presence of a peculiar, red, luminous, flattened stellar halo that may actually correspond to 
tidal debris that in projection appears located closely below the disk plane and that may 
have originated the faint residual light emission reported in these previous studies (as also 
suggested by Zheng et al. 1999). 

What can we learn from the theoretical study of tidal streams in galaxies beyond the 
Local Group? Firstly it must be remarked that theoretical models are severely hampered 
by the impossibility of measuring the kinematics of individual stars at distances of several 
Mpc with present instrumentation. In the case of NGC 5907, the large degeneracies in the 
orbital and structural parameters of the stream progenitor force us to consider our N-body 
model as illustrative only. However, the model is certainly helpful to interpret most of the 
observed features of the NGC 5907 tidal stream and proves that such a complex stream-like 
structure can be explained by the accretion of a single satellite galaxy. 

Despite the degeneracy introduced by the observational limitations, an advantage of the 
panoramic perspective in external galaxies is that it simplifies the detection of stream pieces 
that were lost several Gigayears ago, which strongly increases the number and strength of the 
constraints that can be derived from these systems. For example, (i) the variation of width 
and surface brightness along tidal streams constrains the number of dark matter sub-halos in 
the Milky Way (Ibata et al. 2002; Johnston, Choi, & Guhathakurta 2002; Penarrubia et al. 
2006; Siegal-Gaskins & Valluri 2007). These authors show that the presence of dark matter 
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sub-halos in spiral galaxies would result in a progressive heating of tidal streams as a result 
of close encounters. As Penarrubia et al. (2006) point out, the averaged number of dark 
matter substructures (and, thus, the likelihood of encounters) in a Milky Way-like galaxy 
decreases monotonically since z ~ 2 to the present, which clearly implies that "old" stream 
pieces are more likely to reveal perturbations than recently stripped ones, (ii) Constraints 
on the flattening of the host galaxy potential become significantly stronger as the number of 
stream wraps increase (e.g., Johnston, Sackett, & Bullock 2001; Law, Johnston, & Majewski 
2005; Penarrubia et al. 2005). (iii) Finally, detecting old stream pieces will allow us to study 
possible metallicity gradients within the progenitor galaxy, as it is for example in the case 
of the Sgr tidal stream (e.g., Bellazzini et al. 2006; Chou et al. 2007). 

The detection of an old, complex tidal stream in NGC5907 (14 Mpc away) with modest 
instruments suggests the viability of carrying out extragalactic surveys of substructures in 
the surroundings of spiral galaxies in the Local Volume. It is a relevant question to wonder 
how many stellar streams should we expect in the Local Volume. Cosmological N-body 
simulations predict that the incidence of large satellite galaxy accretion occurring in the 
current epoch is commonly one per spiral galaxy for its inner galactocentric regions, ^ 
50 kpc (e.g., Bullock & Johnston 2005). The presence of such a complex stream at 14 
Mpc from us, easily detectable with present observational capabilities, clearly suggests that 
deep photometric images may reveal a large number of substructures in the stellar halos of 
spiral galaxies. The different stellar streams detected in the Milky Way (Yanny et al. 2003; 
Majewski et al. 2003; Juric et al. 2005; Belokurov et al. 2006a,b) and in M31 (Ibata et al. 
2007) by large photometric surveys supports this expectation. 

Table 2 gives the confirmed stellar stream features discovered up to date. Interestingly, 
the first outstanding characteristic of the tidal streams discovered in external galaxies to 
the present day is there are no detected tidal streams in the Local Group with surface 
brightnesses as high as seen in NGC5907. 

Presently, observational limitations impede the resolution of stars in external galaxies, 
which poses a detection threshold of around fi R ~ 29 mag/ in unresolved light. In contrast, 
star counts can reach a detection level as low as ~33 mag/ □" (see e.g., Belokurov et al. 2006a). 
This explains why we cannot find the counter-parts of the stellar streams detected in the 
Local Group in external galaxies, but it is indeed puzzling that no tidal stream in the Milky 
Way or M31 is remotely as bright as the one detected here. 

In addition, Table 2 also reveals that, with the single exception of the Sgr dwarf, the 
progenitor systems of most of the stellar streams known to the present day remain undetected. 
The most logical explanation is that the progenitor galaxies have been nearly destroyed by 
now. This scenario is supported by the N-body simulations of Penarrubia et al. (2008b), 
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who show that dwarf galaxies undergoing tidal mass stripping have to lose more than 90- 
99% of their initial dark matter halo mass before a stellar tidal stream may start forming. 
Once the host's tidal field begins to strip stellar material, the central surface brightness (E ) 
and the total luminosity (L) of the dwarf galaxy suffers a sharp, monotonic decline so that 
S /S (t = 0) oc [L/L(t = 0)] ab for L/L[t = 0] < 1. Therefore, tidal stream forms the 
likelihood to detect its progenitor decreases until, eventually, the present surface brightness 
falls beyond the detection limits of our instruments. A second possibility that we explore 
in our illustrative model (see Fig. 4) is that the progenitor systems might be hidden behind 
the host disk or bulge. Although this scenario may be applied to some of the systems listed 
in Table 2, it is fairly unlikely that it explains the absence of progenitor galaxies for all of 
them. 

The study of tidal streams in external galaxies is a relatively young field and provides 
direct way of addressing some opened questions on galactic formation and evolution. For ex- 
ample, studying stellar population gradients along tidal streams via deep Hubble Space Tele- 
scope photometry data (see for example Aloisi et al. 2005) will render important constraints 
on the effect of tides on the stellar formation history of dwarf galaxies. The panoramic view 
of tidal streams in external galaxies also offers an excellent opportunity to demonstrate tidal 
stripping of globular clusters formed in satellite galaxies, which may correspond to an impor- 
tant fraction of the globular cluster population of the host, as earlier proposed by Searle & 
Zinn (1978). Searching for associated globular clusters and planetary nebula may also offer a 
means by which kinematics in distant streams could be derived (i.e. from intrinsically bright 
absorption line sources or sources with emission lines), so that kinematical constraints might 
be placed on theoretical models that aim to reproduce tidal streams in external galaxies. 
These models, in turn, may provide reasonable estimates of the flattening of the dark matter 
halos that surround Local Volume galaxies. Ultimately, the ideal scenario would require 
resolving stellar populations at a distance of 10-20 Mpc, which will be feasible in the next 
one or two decades with large ground-based telescopes like the Thirty-Meter-Telescope, or 
spatial missions like the James Webb Space Telescope. 
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Table 1. JOURNAL OF OBSERVATIONS 



uate 


r liter 


Total Exposure Time (s) 


2006 June 6 


Blue 


6480 


2006 June 6 


Clear Luminance 


9000 


2006 June 6 


Green 


1080 


2006 June 6 


Red 


5400 


2006 June 7 


Clear Luminance 


3600 


2006 June 7 


Green 


2160 


2006 June 9 


Clear Luminance 


5400 


2006 August 10 


Blue 


2160 


2006 August 10 


Clear Luminance 


2700 


2006 August 10 


Green 


1080 


2006 August 10 


Red 


1800 



Table 2: Stellar tidal streams detected in nearby spiral galaxies (D < 15 Mpc) 



Name 


Distance (Mpc) 


/i (band) 


Progenitor 


Refs. 


Local Group 


Sagittarius (MW) 




29.6 (V) 


Sagittarius dwarf 


(1) 


Monoceros (MW) 




34.9 (V) 


Canis Major dwarf (?) 


(2),(3) 


Orphan (MW) 




32.4 (R) 


unknown 


(4) 


M31 giant stream (M31) 


0.78 


30.0 (V) 


unknown 


(5),(6) 


External galaxies 










NGC 5236 (M83) 


5.2 


26.5 (R) 


unknown 


(7),(8) 


NGC 5907 


11.2 


26.8 (R) 


unknown 


(9),(10) 


NGC 4013 


12.0 


27.0 (R) 


unknown 


(11) 


NGC 4594 (M104) 


12.4 




unknown 


(8) 


NGC 3310 


14.4 


25.6 (V) 


unknown 


(14) 



References. — (1) Martmez-Delgado et al. (2004); (2) Belokurov et al. (2006); (3) Martin et al. (2004); 
(4) Belokurov et al. (2007); (5) Ibata et al. (2001);(6) McConnachie et al. (2003); (7) dc Jong et al. 
(2207);(8) Malin & Hadley (1997); (9) Shang et al. (1998); (10) this manuscript; (11) Martmez-Delgado et 
al. (2008); (12) Whener & Gallagher (2005) 



Fig. 1. — Image of NGC5907 obtained with the BBRO 0.5-meter telescope. The total 
exposure time of this image is 11.35 hours, co- adding all images obtained in this project 
(see Table 1). The image has dimensions of 18.2 x 27.7 arcmins, which, at the distance of 
NGC5907 is ~ 75 x 115 kpc. For a better comparison with the N-body simulations given 
in Fig. 4, this image is shown east up and north to the right. The linear diagonal feature in 
the lower left corner of the image is spurious light from deflection of a bright star off of the 
edge of the CCD chip. The faint halos surrounding the field stars are due to the pass of the 
Gaussian filter mentioned in Sec. 2.2. 




Fig. 2. — Identified photometric features detected in the image of NGC 5907 obtained with 
the BBRO 20-inch telescope and discussed in Sec. 3.1. For reference, a colour image obtained 
with the same telescope (see Table 1) has been superimposed on the saturated disk region of 
the galaxy. The position of the two different strips used to measure the width of the stream 
(see Sec. 3.2) are also indicated by red rectangles. The yellow lines mark the position of the 
nominal cut- off radius (see Sec 2.2). For comparison purposes, north is to the right and east 
is up. 




Fig. 3. — Gaussian fits to the surface brightness distributions of the stream in two different 
strips (A and B) illustrated in Fig. 2. The FWHM of these Gaussian fits are estimated to 
be ~103 and ~164 arcsec respectively. Note that 1 arcsec at the distance of NGC5907 
corresponds to ~0.068 kpc. 
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Fig. 4. — Illustrative N-body model of the stellar stream detected in NGC 5907. The satellite 
is realized as a King model with an initial mass, King core and tidal radii of M = 2 x 1O 8 M , 
r c = 0.39 kpc and r t = 2.7 kpc, respectively (see Sec. 4.2 for details). Different colours denote 
particles that became unbound after different peri-center passages, whereas black particles 
are those that remain bound to the main system. Our color code is so that cyan, green red 
and blue colours are particles lost during the first, second, third and four orbital periods, 
respectively. For this particular model the orbital period is T r = 0.9 Gyr. 



